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Abstract: We construct charged black hole solutions with hyperscaling violation in the infrared(IR) region in 
Einstein-Maxwell-Dilaton-Axion theory and investigate the temperature behavior of the ratio of holographic shear 
viscosity to the entropy density. When translational symmetry breaking is relevant in the IR, the power law of the 
ratio is testified numerically at low temperature T, namely, n/s ~ T^, where the values of exponent & coincide with 
the analytical results. We also find that the exponent & is not affected by irrelevant current, but is reduced by the 
relevant current. 
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1 Introduction 


In last decade the development of holographic techniques has given new insights into the hydrodynamics with 
strong couplings. One celebrated achievement is the Kovtun-Son-Starinets (KSS) bound for the ratio of shear viscosity 
to entropy density 
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which is considered as a fundamental bound for near perfect fluid with strong interactions. However, in recent 
years counter-examples which violate the KSS bound have been found in holographic literature, including the 
higher derivative gravity [2] [3], anisotropic system [H7] as well as isotropic system without translational invariance 
[S15]. In the latter case, shear viscosity loses its hydrodynamical interpretation because of the non-conservation of 
momentum and is usually defined by Kubo Formula 
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n= lim g mC Teute (w,k=0), (2) 
where C is the retarded Green function of the energy-momentum tensor operator T™ in the dual boundary 


theory. Nevertheless, it is still quite instructive to investigate the temperature behavior of the ratio of shear viscosity 
to entropy density in general holographic models without translational invariance. 

Historically, the breaking of translational invariance is introduced in holography to study the transportation of 
the dual system with momentum dissipation [S] [161/30]. In this setup, new geometries in the IR may emerge and 
often accompany with new scaling relations, leading to new scaling behavior of thermodynamic quantities or Green 
functions with temperature T or frequency w [31139]. 

In particular, by virtue of recent progress in [9] people have learned that when the translation symmetry breaking 
is relevant in the far IR, the ratio exhibits a power law behavior with the temperature 
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which reflects the scaling symmetry emerging in the IR. Moreover, a new bound for the exponent & was proposed as 
k <2 there, which might be supported by a heuristic argument based on the bound for the rate of entropy production, 
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Other holographic models are investigated in [8] [101/12] [14] [15], which also satisfy the bound « < 2. But soon later 
in [13], we find & > 2 is possible when the scaling of Lifshitz [40]42] or hyperscaling violation [22] 43][57] emerges in 
the IR. While the bound for the rate of entropy production is not violated, since the origin of the power law 
should be understood as the nontrivial anomalous dimensions of (2) under the rescaling of the IR solution [9] [13]. 
Furthermore, in [13] we have analytically derived a formula for « in a general Einstein-Maxwell-Dilaton-Axion 
(EMD-Axion) model with spatial dimension d, dynamical critical exponent z as well as hyperscaling violating expo- 


nent 0. Specifically, we find 
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where e? is defined as the ratio of Maxwell term and one of the lattice terms in the Lagrangian. However, in [13] 
only for the case of e? = 0 has this formula been justified by numerical calculation on neutral background. In 
this paper we intend to continuously testify the validity of on charged background within EMD-Axion models. 
Schematically, the case of e? Z 0 can be realized by relevant currents. We will numerically construct specific charged 
backgrounds with ultraviolet(UV) completion and then compute the power law of the ratio of shear viscosity to 
entropy density. As a result, we will show that the formula for the ratio 7/s in which was previously proposed in 
based on the scaling analysis can be justified even for e? Z 0 indeed. 

In this paper, we adopt the statement about ‘(marginally) relevant’ and ‘irrelevant’ in [32]: the current or the 
axion being (marginally) relevant means that the Maxwell or the axion terms is the same order as the curvature term 
and the dilaton potential in Lagrange in the power of the radial coordinate; the current or the axion being irrelevant 
means that the Maxwell or the axion terms is subleading comparing to the curvature term and the dilaton potential. 
Roughly speaking, a field being (marginally) relevant or irrelevant depends on whether or not it would deform the 
IR geometry strongly. 


2 EMD-Axion model and hyperscaling violating metric 


In this section we will present our holographic setup and then outline the logic leading to the formula for the 
exponent & in based on the scaling analysis. The action of a general EMD-Axion model in d+2 space time reads 
as 
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where x;(i = 1,2,---,d) are axions and J(¢),Z(¢),V(@¢) are coupling functions or potential of the dilaton field 9. 
Given above action, the equations of motion can be derived as follows 
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For simplicity, we only consider isotropic solutions with following ansatz 
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where translational invariance is broken by the axions y; but the metric and energy-momentum tensor remain to be 
homogeneous. 


As stated above, we are interested in solutions interpolating between AdS in the UV and hyperscaling violating 
solution in the IR. It can be realized by the process of UV completion |46|]. Firstly, we construct a hyperscaling 
violation solution with running dilaton. Secondly, we modify the local behaviors of potentials to graft the solution 
onto AdS in the UV. Generally, solutions with AdS exist when the dilaton reaches the extremal point of its potential 
and the Lorentz symmetry is maintained [56]. 

Let us find hyperscaling violating solutions firstly. When the potentials behave as 


V(¢)~Voe®%?, J($)ee??, Z(b)~e%?, when $— +00, (9) 


it is found in that scaling solutions with hyperscaling violation exist, whose metric and matter fields read as 
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ds? =r | pue T r2 + m , A-Qr* dt, Car: Xi = kti, (10) 
where k characterizes the scale of breaking of translational invariance and ¢ is called conduction exponent [57]. 
Parameters {z,0,¢,¢,k,Q} in ansatz are determined by equations of motion (7). The solutions are found to be 
classified into four classes [32], depending on the parameters (o, B,"y, Vo) in potentials (9). Among them, the explicit 
solutions with (marginally) relevant axion are shown in Appendix[A] 'The metric in can be deformed into a black 
hole solution 
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From Maxwell equation (7c), we have conserved charge density 
p-v-gZ(o)F". (14) 


It can be shown that there exists a scaling relation 
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We will come back to the UV completion and modify the potentials in > section. 

Now we turn to the study of shear viscosity 7. From Kubo formula (2), 7 can be derived by perturbing (ôg)”, 
9" 0Gxy =h(r)e~*, which is the dual field of operator T". Einstein e uw. aJ give rise to the shear NEUE 
equation 


=, d,(/—gg""O-h(r)) + (g"'w* —m(r)*)h(r) =0, (16) 
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h(r) is required to be regular at the horizon and equal to 1 at the conformal boundary rg. The varying mass is 
supposed to satisfy the condition m(r)? > 0 in the models considered so far. Here, we have m(r)? = J(¢)k?g?*. 
n/s can be obtained by the weaker horizon formula [9] [58] 
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where ho(r) is the solution at w — 0 and r, is the location of horizon. 


Following the analysis presented in [13], one can calculate the exponent & of ? ~ T". Here we present a simpler 
derivation with the use of formula (18). If the axion is (marginally) relevant, by using Einstein equations and 
black hole metric (tip, we have 
2õo(z— 1) 
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where 
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e(r)? is just the ratio of the Maxwell term to one of the axion terms in Lagrangian (6). It goes to a nonzero constant 
in the far IR if the current is also (marginally) relevant, otherwise it goes to zero. Thus at leading order, we just set 
e(r)? =e. If the axion is irrelevant, m(r)? goes to zero in the far IR, then we set M? — 0, which is valid at leading 
order. 


By using and metric (11). we rewrite perturbation equation at w=0 
O,(r! *?» f (r)8.ho(r)) - M? L?r 9?» ho(r) ^ 0. (21) 


Solving this equation we can separately obtain the asymptotical expansion of ho(r) near the boundary and its value 
on the horizon as 
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We have abbreviated the series (=) (n > 1) following (=) and the series (=) (m > 1) 
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Coefficient c should be determined by the boundary condition ho(ra) = 1. r, is the boundary of the region where 
the black hole with hyperscaling violation can be described by (11). When consider the second-order variation 
of the action in (6) over a fixed background with the perturbed metric ho(r), the boundary part of the variation 
with the branch of r?o-?? ( rêt) in is divergent (finite), thus the branch of r?»-?? ( r?7) is non-normalizable 
(normalizable). Then ho(r;)z c (= 

Since the UV completion is taken into account, one should be cautious that r; is just an intermediate scale but 
not the conformal boundary rg. The region between r; and ra is AdS deformed by matter fields. The boundary 
condition requires ho(ra) = 1. As explained in [9], when going from ra to r;, ho(r) decreases monotonously from 1 
to a value T when m(r)? > 0. We have introduced a ‘tunneling rate’ I to characterize how ho(r) tunnels from ra 
to r;. The tunneling rate [ and intermediate scale r; should be independent of T when the scale of T is much less 
than other scales, such as the source of dilaton and k in axion. It is because that scale T is not the dominating scale 
driving the renormalization group(RG) flow from AdS to hyperscaling violating solution. Scale T becomes important 
only when we go into the far IR. Then we have 
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By working out c, we can determine the horizon value 
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By virtue of formula (18), we finally obtain 
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where and have been used. One can also employ UV-IR matching [38] to reproduce the same result 
as what is performed in [I3]. As one can see, when axion is absent or irrelevant, we have m(r)? — 0 at leading order, 
then n/s-- T? [61] [62]. Next we will numerically testify this formula when e? =0 or e? £0. 


3 UV completion and numerical results 


In this section we specifically construct the background interpolating hyperscaling violating solution in the 
IR and AdS solution in the UV with finite temperature and charge density. On one hand, as explained in Appendix 
[A] or in [32], solution can be constructed by choosing exponential potentials (9) with running dilaton. On the 
other hand, AdS can be constructed by finding extremal points of constant dilaton with lorentz symmetry. Here we 
adopt dilaton ¢ to interpolate the UV and the IR solutions. It requires some special settings for the potential V (4). 


3.1 UV completion 


In this paper, we focus only on 0 « d since in this region the entanglement entropy obeys area to volume law, 
which is considered as normal behavior of QFT [13] [53]. In this situation the location of IR is r+ +00. Following 
the discussion in [13] [56], the constraints about (2,0) are reduced to 


(Psonz» Dv (oceani $41). (277) 
Then it leads to ĝo 7 0. We have excluded the two cases of 0 =d and z=  -- 1 which can not be reached by running 
dilaton, as shown in Appendix or [32]. We choose the branch of 6 > 0. From the requirement of the potentials 
(9). one can see that our solution can flow to hyperscaling violating solution in the IR if $ — --oo(r — +00). It 
requires e > 0 in (10). We conduct the UV completion by modifying the potential V(¢) but fixing the other two 
coupling potentials 

J($)-e^*, Z(g)-e*. (28) 
From the analysis in Appendix (A). we find a universal behavior of V(¢) ~ r- in the coordinate of ansatz (19). 
So, when approaching the UV (r — 0), the qualitative behavior of V(¢) depends on the sign of 0. On the side of the 
UV, AdS is allowable if axion and gauge field are turned off and ¢= $, is an extremal point of V(9), where V(¢,) 
stands for cosmological constant. Without loss of generality, we choose ¢, — 0. Then a realistic strategy is modifying 
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In Appendix [A] we have ae = —2?. Then a has the opposite sign of 0. So, as one can see, each V(ó) approaches 
Voe*? when $ — --oo and approaches d(d+1) when $ — 0. Without loss of generality, we have chosen the AdS radius 
to be 1. However, the intermediate behavior is not very important. 

AdS4,» vacuum is always allowable. We choose the first type quantization, and the scaling dimensions about 
the dual source of dilaton ¢ and operator O; are determined by the small ¢ expansion of V(ó). When 6 = 0, as 
dilaton is massless, we have A,(o; — 0 and Ao, =d+1, which is marginal deformation. We expect it to be marginally 
relevant to drive the solution away from AdS, as ¢=0 is not a stable point when axion and gauge field are turned 
on. When 6 Z 0, as V(¢) — d(d-- 1) - 29? 4----, we have Apo) —1 and Ao, — d, which is relevant deformation. 


3.2 Numerical calculation and results 


We use the following ansatz for numerical calculation 


2 1 —S(u) 442 du? : 2 
d? =- (-a-uvw dt +e), 
ó—ó(u) xyi=ka;,, A-(1—u)A(u)dt. (30) 


'The conformal boundary is located at u — 0 while the horizon at u — 1. The temperature and the entropy density 
are T = -U(1)e 50/7 and s— 4m. 

The AdSy,2 vacuum corresponds to U = 1, S = ọ = A = 0. Boundary conditions at the horizon are regular 
conditions. Boundary conditions at conformal boundary should satisfy the scaling dimensions, which depend on the 
potential V(¢) as well as the value of 0. 


Explicitly, the asymptotic expansions near the conformal boundary are 
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where p is the chemical potential and e79 has been set to 1 by rescaling t. The different boundary conditions of 
ol(u) come from the different choices of V(¢) in (29). We can work on either grand canonical ensemble or canonical 
ensemble. 


3.2.1 Grand canonical ensemble 


In grand canonical ensemble, we control the value of chemical potential ju. 

When 0 = 0, the boundary conditions at conformal boundary are U(0) — 1, S(0) =0, A(0) =p, ¢(0) =A. We 
work in the unit of k. The dimensionless quantities parameterizing the family of black hole solutions are is KAL 
When 0 #0, the boundary conditions are U(0)=1, S(0)=0, A(0)= u, ¢'(0)= A. The dimensionless quantities are 

T H X 

, We T construct the interpolating solutions for ¢ > 0. When dropping down T/k, we should fix the values 
of (4,A) (for 0 20) or (£,2) (for 0 Z 0) within an appropriate region respectively, in order to reach hyperscaling 
violating solution in the IR at low T/k. 

The dimensionless entropy density and charge density are s/k^ and p/k*. We calculate n/s by using and 
found 2 < i for all the time, because of the breaking of translational invariance. 

At high T/k, the scaling relation is controlled by AdS in the UV, which gives rise to the power laws of s ~ T7 and 
€ T°. On the other hand, at low T/k, the hyperscaling violation emerges in the IR. The power law of s ~ ge 
and ? e T" are observed in numerical results. 

It is worthwhile to notice that p/k^ converges to a nonzero constant at low T/k. When approaching the IR, 
e*(u) converges to a nonzero constant for class I but goes to zero as some power of radial coordinate u for class II. 
Similarly, at low T/k, the horizon value e? = e?(1) converges to a nonzero constant for class I but goes to zero as 
some power law for class II, whose exponent is shown in Appendix [A] For the same (d,z,0), it is observed that the 
appearance of a nonzero e? always reduces the exponent « of 7 ~ T^, which is consistent with the property that & 
decreases monotonously with e? in when & 0. 

We conduct the numerical calculation for d = 2 and 0 = si 0, —4 as representives of three cases, namely 0 < 0, 0=0 
and 0 « 0 « d. Different values of are chosen to represent class I or class II. The specific results are shown in Figure 


All the numerical results match the analytical ones. 
3.2.2 Canonical ensemble 


In canonical ensemble, we control the value of charge density p. The t component of Maxwell equations can 
be replaced by (14). 

When 0 = 0, the three-parameter family of solutions is characterized by { T, zr, Aj. When 8 Æ 0, the one 
is characterized by (7, 4, >}. When lowering down T/k, we fix the other two dimensionless parameters within 
an appropriate region. Then u/k converges to a nonzero constant at low T/k instead of p/k? in grand canonical 
ensemble. The behaviors of 7/s and e? are similar to those in grand canonical ensemble. By using the method in 
Appendix [A] we can foresee the value to which e? converges in the IR at low temperature for class I. 

In a parallel manner, we conduct the numerical calculation for d = 3 and 0 = —6, 0, 2 to represent the three 
regions of 0 < 0, 6 — 0 and 0 « 0 « d. The specific results are shown in Figure [4] [5] and [6] All the numerical results 
match the analytical ones. 


4 Conclusion and Outlooks 


In this paper we have numerically constructed charged solutions with emerging hyperscaling violation in EMD- 
Axion model and investigated the temperature behavior of the ratio of shear viscosity to entropy density. We have 
found that the relevant axion, which breaks the translational invariance, leads to the power law of 7/s ~ T^. In 
particular, the relevant current reduces the exponent « indeed. This reduction is characterized by the quantity e?, 
which can be derived from dimensionless conserved charge density p/k^. While, irrelevant current does not affect the 
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^y — 0.6 belongs to class I; y = 2 belongs to class II. Solid lines denote numerical results and dashed lines denote 
analytical results from and (26). Upper-left: Plot of e?(u) at rather low temperatures, where it converges 
to a nonzero constant in the IR for y=0.6. Upper-right: Log-log plot of p/k? as a function of T/k, where p/k? 
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exponent & since e? > 0 in the far IR at low temperature. Our analytical results for the exponent & coincide with our 
numerical calculation, indicating that our proposed formula for « in is robust at least for generic backgrounds 
within EMD-Axion model. 

Especially, our results for Lifshitz case verify that & > 2 can happen in d > 2, indicating that hyperscaling violation 
is not the essential ingredient leading to the exponent & 7 2. Moreover, as conjectured in [13], the upper bounds for 
k coincide with the behaviors of entanglement entropy. 

Analytically, it is possible that irrelevant current or axion can affect the temperature behavior of 7/s at subleading 
orders. One should consider their backreaction to the background, then solve the shear perturbation equation (16). 
It is related to the issue whether temperature T is still the unique scale in entropy production, as it is when axion 
is relevant. It is worthy of investigation in the future. 

We are very grateful to Blaise Goutéraux, Peng Liu and Xiangrong Zheng for helpful discussions and correspon- 
dence. 


A The classification of IR solutions 


We focus on relevant axion solutions, otherwise 7/s just converges to a non-zero constant at low temperature at 
leading order. We require the solutions should have positive specific heat, and the temperature deformation is the 
only allowed relevant deformation. In addition, we give an extra requirement of 0 « d. The scaling solutions have 
been obtained and classified in [31] [32]. 


A.1 Class I: marginally relevant current, marginally relevant axion 


If both the Maxwell and the axion terms are the same order in the power of the radial coordinate as the curvature 
term and the dilaton potential V($) in Lagrange (6), scaling solutions obtained form a one-parameter family 
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which can be parameterized by k in coordinate of (10). The charge related quantities are 
ó2Q? — k?(0—dz)-4-2Vo(z—1) 
2 _ 0% .— 
xls. dx AI , (33) 
à = wa _ k?(0—dz)+2Vo(z-1) (34) 
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The mode analysis in [32] indicates that there are three pairs of conjugate modes summing to 69. Two pairs 
are degenerate with £,.. = 6, = 0 and 81,4 = B3. = ĝo. f, rescales the time. D, is temperature deformation 
and responsible for creating a small black hole (11). B», changes k and shifts the solution along the one-parameter 
family. 82,4 changes the chemical potential and belongs to the transformation of gauge symmetry. The expression 
for the last pair of 63.4, is too tedious to show here. We require that 81, is relevant and 83 is irrelevant, with IR 
located at r — oo. Then the final allowed parameter space here is found to be p? > 0 and (27). 


M -attgrr gt, Z(o)F"* 
(85 x)* 
we can use it to connect UV with IR and determine the solution in the one-parameter family by using (33). Finally, 


€? can be obtained from £ by using (34). which is rather convenient in canonical ensemble. 

At zero temperature, one can integrate the three modes of £5.., 85. and (3, to the UV and adjust them to 
satisfy the boundary conditions specified by (A, p, u} at the conformal boundary. A finite temperature solution can 
be driven by £,.,. 


Since the quantity 4 = is conserved and invariant under coordinate transformation within B. 


A.2 Class II: irrelevant current, marginally relevant axion 
If only the Maxwell term turns to be subleading, a single scaling solution at leading order is obtained as 


Fa, wee pee 
| óo(dz — 0) a. 2Va(z— 1) 
EN md (35) 


L? 
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There are three pairs of modes, in which two pairs sum to do. The first pair is 64,- = 0 and £,, = ĝo which are 
rescaling of time and temperature deformation. The second pair is relevant £5, and irrelevant 62. The third pair 


is gauge field modes with 


20 
A(r) = Ai - Ayr* 7, C =d-0+ — — ey. (36) 


d 
¢-—d+0 


Mode of A; is irrelevant when óo(C4-d—0) <0. We find e?(r)  r*74^** and p~T°,e2 ~T -* , when T > 0. Similar 
to class I, we can integrate B5... and gauge field modes to the UV. 
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